BUSEMANN FUNCTIONS AND EQUILIBRIUM MEASURES 
IN LAST PASSAGE PERCOLATION MODELS 



ERIC CATOR AND LEANDRO P. R. PIMENTEL 



Abstract. The interplay between two-dimensional percolation growth models and one-dimensional 
particle processes has been a fruitful source of interesting mathematical phenomena. In this paper 
we develop a connection between the construction of Busemann functions in the Hammersley last- 
passage percolation model with i.i.d. random weights, and the existence, ergodicity and uniqueness 
of equilibrium (or time-invariant) measures for the related (multi-class) interacting fluid system. 
As we shall see, in the classical Hammersley model, where each point has weight one, this approach 
brings a new and rather geometrical solution of the longest increasing subsequence problem, as well 
as a central limit theorem for the Busemann function. 



1. Introduction 

In the middle of the fifties H. Busemann |6j introduced a collection of functions to study geomet- 
rical aspects of metric spaces. These functions are induced by a metric d, and by a collection of rays 
(semi-infinite geodesies) as follows: the Busemann function b m (-), with respect to a ray (ro(r)) r >o, 
is the limit of d(w(r),w(Q)) — d(zu(r), •) as r goes to infinity. Along a ray w the metric d becomes 
additive. By using the triangle inequality, this implies that the defining sequence is nondecreasing 
and bounded from above, and so it always converges. Using analogous considerations, one can 
construct Busemann functions over spaces equipped with a super-additve "metric" L (one needs 
the reversed triangle inequality) . In this work we are particularly interested in geometrical aspects 
of the following stochastic two-dimensional last passage (super-additive) percolation model: let 
P C R 2 be a two-dimensional Poisson random set of intensity one. On each point p 6 P we put 
a random positive weight uj p and we assume that {w p : p € P} is a collection of i.i.d. random 
variables, distributed according to a distribution function F, which are also independent of P. 
When F is the Dirac distribution concentrated on 1 (each point has weight 1; we will denote this 
F by 5i), then we refer to this model as the classical Hammersley model (Aldous & Diaconis [I]). 
For each p,q 6 IR 2 , with p < q (inequality in each coordinate, p 7^ q), let II(p,q) denote the set 
of all increasing (or up-right) paths, consisting of points in P, from p to q, where we exclude the 
starting point p. In this probabilistic model, the "metric" (or last-passage time) L between p < q 
is defined by 

L(p,q) := max { V o; p /} . 
roeil(p,q) *—? 

Then L is super-additive, 

£(p,q) > £(p,z) + L(z,q) . 
When we consider a path w from p to q consisting of increasing points (pi, . . . , p n ), we will view 
w as the lowest increasing continuous path connecting all the points, starting at p and ending at q, 
and then excluding p. This way we can talk about crossings with other paths or with lines. A finite 
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geodesic between p and q is given by the lowest path that attains the maximum in the definition 
of L(p, q), which we will denote by w(p, q) (this is well defined for any ordered pair (p, q), even if 
we do not specify the order). 

In [9], using methods developed by Newman and co-authors in [IT] and [2] and also applied to 
the classical Hammersley process by Wiithrich in |20| . it is shown that these finite geodesies can 
be extended to semi-infinite a-rays by moving one endpoint appropriately to infinity. An a-ray 
starting at p 6 S 2 , denoted by zu a (p), is a semi- infinite geodesic that starts at p and moves to 
infinity in the direction a := (cos a, sin a), for a £ (7r,37r/2). It turns out that for fixed a, with 
probability 1, each p £ M? is the starting point of a unique a-ray, and two a-rays will always 
coalesce eventually. In Section [2] we will state precisely the theorems for existence and coalescence 
of a-rays. An important tool in the development of these concepts is the following result, known 
as the shape theorem: there exists a constant 7 = j(F) such that for all x, t > 



In Section [3] we will use the a-rays to construct the Busemann function B a : if x,y G M 2 , let 
c = c(a,x, y) be the coalescence point of the two a-rays starting at x and y. Then 



It is not hard to see that the distribution of B a is invariant under translation and that B a is 
additive: B a (x,z) = B a (x,y) + B a (y,z). Its most important property, however, is the following 
connection with an associated interacting particle process: if we define the measures v l a on R by 



then the family of random measures \v l a : t £ R} forms a Markov process, and its evolution 
corresponds to the evolution of an interacting particle process (maybe interacting fluid process is a 
better name) which is a natural extension of the classical Hammersley interacting particle process, 
where F = S±. This connection turns out to be the key idea of this paper. On the one hand, it 
allows us to prove existence, uniqueness and ergodicity (mixing) of the equilibrium measures of the 
interacting fluid process. Furthermore, it gives us a natural way to prove a strong law of large num- 
bers for the second class particle and to define a multi-class fluid system with a countable number 
of classes, which is a new result even in the classical Hammersley process. Also, the classical result 
that 7(^1) = 2 ( (jl.ip for the classical Hammersley process) follows easily from our methods. On 
the other hand, it implies a central limit theorem for the Busemann function in the classical model. 
This result also shows a phase transition from a Gaussian limit distribution, on the square-root 
scale, to (zero-mean) Tracy- Widom type limit distribution, on the cube-root scale, at the critical 
angle a — it. This transition from square-root to cube-root scaling for Busemann functions was first 
conjectured in 2001 by Howard & Newman p3] (see also Wiitrich |20|). 

In Section [8] we describe how to define a Busemann function in the last passage percolation 
on the lattice 1? with iid weights on the lattice points. It is an important open question for a 
long time already, how to prove in this general setup that the shape function, see (jl.ip . is strictly 
curved, a fact we need to define our a-rays. However, if we restrict ourselves to exponential (or 
geometric) weights, in which case the last passage percolation is an alternative description of the 
totally asymetric exclusion process (TASEP), we know the shape function, and it is indeed strictly 
curved. Therefore, we can define the Busemann function and we find a similar connection to the 



(1.1) 




B a (x,y) = L(c,y 



)-L(c,x). 



u t a ((x,y]) = B a {(x,t),{y,t)) 
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(known) equilibrium measures, which allows us to prove analogous results. 

We feel that the connection between the Busemann functions and the equilibrium measures gives 
us an important new tool to study last passage percolation and the corresponding fluid processes. 
In an upcoming paper, we will use the Busemann function and the results from this paper to de- 
termine the asymptotic speed of a second class particle, given a deterministic, rarefaction initial 
condition, in the classical Hammersley process and in TASEP. Furthermore, we have strong indi- 
cations that the Busemann function can help us establish the cube-root behavior of the length of 
a longest path, and the fluctuations of the longest path. 

Overview. In Section [2] we will state the theorems, that we will use further on, about the limit 
shape and the existence and coalescence of a-rays, and then give precise references for the proofs. 
In Section [3] we define the Busemann function and give its most important properties. In Section 
U]we introduce the Hammersley interacting fluid process and establish the connection between the 
Busemann function and the equilibrium measures. In Section [5] we prove uniqueness and ergodicity 
(mixing property) of the equilibrium measures, and we show local convergence to the equilibrium 
measure in case of a rarefaction fan. In Section [6] we prove the central limit theorem for the 
Busemann function in the classical model. In Section we show how we can define a multi-class 
system with a countable number of classes, and we establish the strong law for a second class 
particle. In Section [8] we state the analogous results for the TASEP. 



A key notion in this paper will be an a-ray: for each angle a 6 (7r,37r/2) and for each point 
x£i 2 , Tt7 a (x) is the lowest continuous down-left path through an ordered sequence (pi)«>o m M 2 , 
with po = x, pi G P and pj > p^- whenever i < j. Furthermore, vj(pj,Pi) C w a (x) (every part of 
the path is a geodesic), and finally we must have that 



A crucial step for the existence of a-rays is the following shape theorem: set = (0, 0), n = (n, n), 



2. Shape function and «-rays 



(2.1) 




Pi 



F(x) = P(uj p < x) and 7 = j(F) = sup 



E(L(0,n)) 



> 0. 



n>l 



n 



Theorem 2.1. Suppose that 



00 



(2.2) 



y 7 ! - F(x) dx < +00 . 



Then 'y(F) < 00 and for all x, t > 0, as r — > oo, 




and 




r 



r 



Further, if (|2.2p is strengthened to: there exists a > such that 



(2.3) 




then there exist constants cq, c\, C2, C3, C4 > such that for all r > cq 
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for u G [ c^^fr log 2 r , C4r 3 / 2 log r ] . 

For proofs see Theorem 1.1 and 1.2 in [9]. Theorem 12.11 shows that L has a curved limiting shape, 
mainly due to the invariance of the Poisson process under volume preserving maps: if x, t, r > 
and pel 2 , then 

(2.4) L(0,{x,t))=L(p,p + (rx,t/r)) . 

This is because under this map, the distribution of the Poisson process does not change, and the 
up-right paths are preserved. The almost sure convergence is a standard consequence of the sub- 
additive ergodic theorem, once we have a bound on E(L(0,n)), linear in n. 

The following theorem gives us existence and coalescence of a-rays. Before we state the theorem, 
we shall define what we mean by convergence of paths: we say that a sequence of paths w n converges 
to w, and denote lim n _ i . 00 w n = w, if for all bounded subsets B C M 2 there exists no such that 
w n n B = w n B for all n > n Q . 

Theorem 2.2. Assume (|2.3p , so for some a > 0, 

POD 

/ exp(ax) dF{x) < +oo . 
Jo 

Then for fixed a G (it, 3tt/2) the following holds with probability one: 

(1) For each x G M. 2 there exists a unique a-ray starting from x, which we denote by w a (x). 

(2) For any sequence (z n ) ra >o of points in M 2 with ||z n || — > oo, 

if lim n = (cos a, sin a) then lim ro(x, z n ) = ro a (x) . 

(3) For all x, y G M 2 there exists c = c(a,x,y) such that w a (x) and w a (y) coalesce at c: 

ct7 a (x) = w(x, c) U w a (c) and w a (y) = w(y, c) U w a (c) . 

The proof of Theorem 12.21 is based on a method introduced by Newman |17| that can be applied 
in a wide percolation context. An outline of the proof is given in [9], Theorem 2.4 and 2.5. A 
detailed proof of the same theorem, but now restricted to the Hammersley classical model, can be 
found in [20]. 

3. The Busemann function 

Using the concept of a-rays, we will study the function £> Q (x,y), which is defined by taking the 
first coalescence point c = c(a,x,y) between the a-ray that starts from x and the one that starts 
from y (remember that these two rays coalesce), and setting 

(3.1) 5 Q (x,y) = L(c,y)-L(c,x). 

Note that if we take a different coalescence point c', then c > c' and they both lie on a geodesic. 
Since L is additive on a geodesic, we get L(c', x) = L(c' , c)+L(c, x), which shows that the definition 
of -B a (x, y) does not depend on the choice of the coalescence point. Let (z ra ) n >i be any unbounded 
decreasing sequence that follows direction (cos a, sin a), and let c(z n , x, y) denote the most up-right 
coalescence point between ro(z n ,x) and ro(z n ,y). By Theorem 12.21 with probability one, there 
exists no > such that 

(3.2) Vn>n c(z n , x, y) = c(a, x, y) and L(z n ,y) - L(z n ,x) = B a (x,y) . 
Therefore, in geometrical terms, B a (x, •) can be seen as the Busemann function along the ray zu a (x). 
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Some properties of B a are summarized in the following proposition. The proofs are relatively 
straightforward, and can be found in the Appendix. 

Proposition 3.1. Define the Busemann function B a as above, for a E (tt, 3tv/2). 

(1) The distribution of the function B a is translation invariant: Vp£l 2 

B a (-+p,- + p) =B a (;-). 

(2) B a is anti- symmetric and additive: V x, y, z £~M? 

B a (x,y) = -B a (y,x) and B a (x, z) = B a (x, y) + B a (y, z). 

(3) For any (x,t) £ R 2 , 

B 57r/2 _ a (0,(x,t))=B a (0, (t,x)). 

(4) If x < y and x/y, then 

B a (x,y)>0 and < K(B a (x, y)) < +oo. 

(5) Fix x,y 6 R 2 and p, q € K 2 such that p, q > 0. The function A t-t B a (x + Ap,y + Aq) is 
cadlag in A £ R. 

As mentioned in the Introduction, the most important aspect of B a is a Markovian structure 
described in the following proposition. 

Proposition 3.2. For all s < t and x £ M we have 

B a {(0, s), (x, t)) = sup {B a {(0, s), (z, s)) + L{(z, s), (x, t))} . 

z<x 

Proof: Without loss of generality we can take s = (and therefore t > 0). Define Z a = Z a (x, t) & M 
as the crossing-point of the a-ray starting at (x, t) with the x-axis. Clearly, Z a < x and 

B a (0,(x,t)) = B ol (0,(Z a ,0))+B a ((Z a ,0),(x,t)) 
= B a (0,{Z a ,0)) + L{(Z a ,0),{x,t)). 

The last equality follows from the fact that (x,t) and (Z a ,0) are lying on an a-ray. This means 
that it is enough to prove that for all z < x, 

(3.3) B a (0, (z, 0)) + L((z, 0), (x, t)) < B a (0, (Z a ,0)) + L((Z a , 0), (x, t)). 
Suppose p is a coalescence point of the a-rays starting at 0, (x,t) and (z, 0). Then 

B a (0, (z, 0)) = L(p, (z, 0)) - L(p, 0) and B a (0, (Z a ,0)) = L(p, (Z a ,0)) - L(p, 0). 

Furthermore, since p, (Z a ,0) and (x,t) are elements of Tr a (x,t), we know that 

L(p,(z,0))+L((z,0),(x,i)) < L(p,(x,*)) 

= L(p,(Z a ,0))+L((Z a ,0),(x,t)). 

From this, (|3.3p easily follows. 

□ 

Define the positive measures vf on M, using Proposition 13.11 by 

(3.4) y]) :=S a ((x,t),(0,0)) - B Q ((y,t), (0,0)) = 5 a ((x, t), (y, t)) V x < y. 

Proposition 13.21 then shows that the process 1 1— > vf is a Markov process: the future evolution of vf 
depends on the Poisson process in the upper-half plane M x (t, oo) and on the present value of ff, 
not on the past of the process (which is of course independent of the Poisson process in M x (t, oo)). 
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v{x) 



Note that the distribution of vf does not depend on t (by Proposition 13 . 1 H fTj) ) . so this distribution is 
an equilibrium (or time invariant) measure for the underlying Markov process. In the next section 
we will describe the generator of this Markov process, which will be an extension of the classical 
Hammersley interacting particle process. 

4. The Hammersley interacting fluid system 

It is well known that the classical Hammersley model, where all weights are 1, described in Aldous 
& Diaconis [I], has a representation as an interacting particle system. The Hammersley process 
with random weights has a similar description, although a better name might be an interacting 
fluid system. We start by restricting the compound Poisson process {u p : p G P} to R x R + . 
Then we choose a positive, locally finite measure v defined on R. Usually, these measures will be 
purely atomic, but this is not necessary. To each measure v we associate a non-decreasing process 
v(j) defined by 

u([0,x\) for x > 
-u((x, 0)) for x < 0. 

Note that v{-) is a cadlag function. Although the details are a bit cumbersome, all the results we 
will show can be extended, mutatis mutandis, to the case where v{x) = — oo for x < 0, which would 
correspond to a non-locally finite measure with an infinite fluid density to the left of 0. This is a 
quite natural starting condition, but we will not use it explicitly in this paper. 

The Hammersley interacting fluid system (M" : t > 0) is a stochastic process with values in 
the space of positive, locally finite measures on R. Its evolution is defined as follows: if there is a 
Poisson point with weight w at a point (xo,t), then M^{{xq}) = M"_{{xq\) + oj, and for x > xq, 

(4.1) M?{(x ,x]) = (AC((x ,x])-w) + . 

Here, M t 1 is the "mass distribution" of the fluid at time t if the Poisson point at (xQ,t) would be 
removed. To the left of xq the measure does not change. In words, the Poisson point at (xo,t) 
moves a total mass u to the left, to the point xq, taking the mass from the first available fluid to 
the right of xq. See Figure Q] below for a visualization, in case of atomic measures, of the process 
inside a space-time box. In this picture, restricted to [0, x], the measure v consists of three atoms 
of weight 5, 3 and 7. The measure M", 2 consists of three atoms of weight 1, 4 and 6, while at time 
t, it consists of one atom with weight 7. 

It is not true that the evolution Mj is well defined for all measures v (e.g. if we start with a finite 
number of particles to the left of 0, every particle would be pulled instantaneously to — oo). In 
this paper we follow the Aldous & Diaconis pQ graphical representation in the last-passage model 
(compare to the result in the classical case, found in their paper): 

Proposition 4.1. Let N be the set of all positive, locally finite measures v such that 

(4.2) liminf^^>0. 



y^r—oo y 



For each v £ N , the process defined by 



(4.3) L u (x,t) :=sup{v(z)+L((z,0),(x,t))} (x £ R, t > 0) 
is well defined and the measure 

(4.4) M?((x,y]) := L u (y,t) - L u (x,t) , 
evolves according to the Hammersley interacting fluid system. 
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Figure 1. Example of the Hammersley interacting fluid process 



Proof: It follows immediately from the definition that L v {x,t) is increasing in x and t, even if L v 
would not be finite everywhere. This implies that if we can prove that L v is finite on for example 
7L x Z + with probability one, then almost surely, L v is finite everywhere. Therefore, we need only 
prove for any fixed point (x,t) that L u (x,t) is finite with probability one. Use Theorem 12.11 and 
(H3D to see that 

EL ((z, 0), 0, t)) ~ 7i/t(s - z) 0^-oo). 

The Markov inequality and (|4.2p then give us the desired result. 

We can show that x h-> L v ((z,0), (x,t)) is cadlag using the fact that for y > x, 

L((z,0),(y,t)) < L((z,0),(x,t))+L((x+,0),(y,t)) and UmL((x+, 0), (y, t)) = 0, 

where x+ means that you are not allowed to use a possible Poisson point directly above x. Therefore, 
Mf is indeed a locally finite measure on R. To see that follows the Hammersley interacting 
fluid dynamics that we have just defined, suppose that there is a Poisson point in (xq, t) with weight 
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lj. For x < xo, this Poisson point has no effect, so M t " = Mf_ on (—00, xo). Clearly, 

M?({xo})=MZ_({x }) + u. 

If x > xo, then the longest path to (x,t) that attains the supremum in (|4.3p can either use the 
weight in (xo,t), which would give L u (x,t) = L v (xo,t) = L u (xo,t—) + u, or it could ignore the 
weight in (xo,t), which would give L u (x,t) = L v (x,t—). This proves that 

M?((x ,x]) = L v {x,t) - L v {x ,t) 

= ma,x(L u (x ,t-) + oj, L u (x, t—)) - L u (x ,t-) - u> 
= (MU(xo,x])-u) + . 

□ 

Assume that we have a probability measure defined on M and consider v G M as a realization 
of this probability measure. We say that v is time invariant for the Hammersley interacting fluid 
process (in law) if 

Af t " = M^ = v for all t > . 

In this case, we also say that the underlying probability measure on J\f is an equilibrium measure. 
Let a E (ir, 3tt/2) and define the measure in M 

v a ((x,y])=B a {(x,0),(y,0)) for x<yeR. 

Compare this to (|3.4p {y a = Uq). It was this interplay between the longest path description and 
the equilibrium interacting particle system that proved very fruitful in the results for the classical 
Hammersley process in Cator & Groeneboom [8J. We will attempt the same in the interacting fluid 
system, but since the equilibrium solution is not explicitly known, we needed to develop new tools 
and ideas, which in fact also had interesting applications for the classical case. Of course, as an 
immediate consequence of Proposition 13.21 we do have the following 

Corollary 4.2. The random measures v a = Vq (j3.4j) are all equilibrium measures for the Ham- 
mersley interacting fluid process. 

5. ERGODICITY AND UNIQUENESS OF THE EQUILIBRIUM MEASURE 

Let a £ (it, 3ir/2) and define the measure in N 

v a ((x,y]) =5 a ((x,0),(y,0)) for x<yeR. 

Compare this to (|3.4p [y a = Vq). To prove the next theorem about time invariance, we need to 
define the following exit points: 

(5.1) Z a (x,t) is the crossing point of zu a (x,t) and R x {0}. 
This was already used in the proof of Proposition 13.21 Analogously, 

(5.2) Z*(x,t) is the crossing point of w a (x,t) and {0} x R. 
Finally, if v E M, we define 

(5.3) Z v (x, t) = sup{z < x : L u (x,t) = v(z) + L((z,0), (x,t))}. 

This says that Z u (x,t) is the right-most point where the supremum in the definition of L u (x,t) 
is attained. This definition is slightly subtle. The fact that the supremum in the definition of 
L u (x,t) is attained, relies on the fact that on a compact interval, the sum of a non-decreasing 
right-continuous function iy{z)) and a non-increasing Zefi-continuous function {z 1— > L((z,0), (x,t))) 
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attains its maximum. Now we have to show that for all (y, s), the supremum over z < y can actually 
be restricted to a compact set. For a countable set of (x, t)'s this can be done as in the proof of 
Proposition using (|4.2p and Theorem 12. II for L. Then we can conclude the desired compactness 
property for all (y, s) by using the inequality 

v{z) + L((z, 0), (y, s)) < u{Z v {x, t)) + L{(Z v {x, t),0), (y, a)) V z < Z u (x, t),y>x,s< t. 

We have the following proposition. 

Proposition 5.1. For all x E K and t > 0, we have 

Z a (x,t) = Z Va (x,t). 

Proof: From (|3.3|) we immediately get that Z Ua (x,t) > Z a (x,t). Now suppose Z Ua (x,t) > Z a (x,t). 
Since L(Z a (x, t), (x, t)) = L(Z Ua (x,t), (x,t)) + B a (Z Ua (x,t), Z a (x,t)), we could have chosen the a- 
ray starting at (x,t) through the point Z Ua (x,t), which would make it strictly lower than w a (x,t), 
going through Z a (x,t). □ 

Since the ct-ray w a {Q) has asymptotic direction (cos(a), sin(a)) with probability 1, we can easily 
see that 

(5.4) lim P(\Z a (t, (tan a)t)\ > et) = lim P(\Z*(t, (tana)t)| > et) = . 

t— >oo t— >oo 

However, we can also control Z v for more general za Compare the following lemma to Lemma 3.3 
in Ferrari, Martin & Pimentel 

Lemma 5.2. Suppose v £ N . Assume that 

(5.5) liminf — ' — — > ^-Vtana and that lim sup — - < ^ V tan a . 

z^-oo —z 2 z^oo Z 2 

Then, with probability one, 

Z v (t, (tana)i) 
hm = . 

t— >oo t 

Proof: Using the transformation (|2.4|) . we can assume, without loss of generality, that a = 5tt/4. 
The proof of this lemma is based on the following estimate for the shape function f(x,t) = r )\fxi: 
for s £ [0, t] 

2 

(5-6) f( t ,t)-f(t-a,t)>la + lj. 

Now, fix e > and suppose that Z u (t) := Z u (t,t) > et + 1. Then there exists k E [et, i] H N and 
n = [tj such that 

+ l)+L((fc,0),(n + l,n + l)) -L((0,0),(n,n)) >0. 

By adding f(n, n) — f(n — k,n) — ^k to both sides of the above inequality and applying (j5.6|) . we 
get 

u(k + l)-^k + 
L((k, 0), (n + 1, n + 1)) - f(n -k,n) + 

f(n,n)-L((0,0),(n,n)) > 1^ 

(5.7) > Je 2 n. 
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By (|5.5p . we know that 

r !<!+Hf 

iim sup sup < (J . 

Therefore, by (15. 7p . if the set {t > : Z v (t,t) > et + 1} is unbounded, it follows that for some 
small r\ > 0, the events 

{|L((fc,0),(n + l,n + l)) -/(n-fc,7i)| > r?n} 

happen infinitely often for n > 1 and < fc < n. Using Borel-Cantelli, this will have zero probability 
if for all rj > 0, 

oo n 

(5.8) S P °)» ( n + " + X )) ~ ~ n )l > < +°°- 

n=l k=0 

Theorem 12.11 gives us some control on the fluctuations of L about its asymptotic shape. Note that 
for n big enough, and < k < n, we have 

P(|L((A;,0),(n + l,n + l)) - f(n-k,n)\ > r/n) = 

= P (|L(0, v / (n + l-A;)(n + l)(l, 1)) - /(n - fe,n)| > r]nj . 

Define r = y/(n + 1 — k)(n + 1) and it = n 2 / 3 . If we choose re large enough, we can make sure that 
u < rjn/2 and 

\f(n — A:, re) — jr\ < rjn/2. 

Also, for n large enough, we have that u G [ 03^^ log 2 r , C4r 3//2 logr] (see Theorem 12. ip . This 
implies, using Theorem l2.ll that there exist c±, C2 > such that for n large enough and < k < n, 

P(|L((fc,0),(n + l,n+l)) -/(n- k,n)\ > rjn) < P(|L(0, (r,r)) — jr\ > uj 

< c\ exp ^— C2re 1 ^ 6 / log(re' 

This clearly proves (|5.8p . The proof that the set {t > : Z u (t,t) < —et — 1} is bounded with 
probability 1 follows the same line. □ 

Now we can prove the most important result of this Section. 
Theorem 5.3. If we start the Hammersley interacting fluid system with v a then 

M v t a = u a for all t > . 
The process x 1— > f a {x) is stationary and ergodic and its intensity is given by 

(5.9) Ei/ a (l) = ^p-Vtaaa. 

Finally, consider a random v £ N , which is time invariant, and which defines a stationary and 
ergodic process on R. Define a S (7r,37r/2) by 

2 , x 2 

a = arctan ( — — 
7 (F) 



Then v = v a . 
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Proof: The first statement is an immediate consequence of Proposition 14.11 and Proposition 13.21 
The fact that x i— > v a (x) is stationary and that < E(f a (l)) < +oo follows directly from Propo- 
sition 13.11 Ergodicity will follow from Proposition 15.61 but won't be needed for the rest of this 
proof. Fix a £ (tt, 37r/2) and set p = p(a) := Vtana. Since the model is invariant under the map 

(x,t) -> (px,t/p), 

Ei/ a (l) = Eu 5n/4 (p) = Ez^ 57r/4 (l)Vtana. 

Now, for all t > 0, 

L 5jr/4 (t,t) = i^ r/4 (z v ^(t,t)) +L((Z„ W4 (M),0),(M)) • 
By Proposition EH Z„ 6 _. /4 (i, i) = Z 57r/4 (i,i) =: ^ 57r / 4 (i). If Z 5n/4 (t) > then 

0<L l/W4 (t,t)-L((0,0),(t,t)) < v 5n/i (Z W4 (t)) . 
On the other hand, if Z 57I y 4 (t) < then Z^,^{t) > 0. Now define 

(5.10) i/*(x) = £ Q (0,(0,x)). 

From the additivity of the Busemann function, we know that 

L 57r/4 (M) = ^ /4 (^ /4 (i,t)) + L((0.Z£ r/4 (i,t)).(M)) . 
Finally we obtain that, denoting Z*_y 4 (i, t) =: Z^,^(t), 

< L U5w/4 (t,t)-L((0,0),(t,t)) < max{^ /4 (Z 5V4 (t)),^ /4 (Z* V4 (t))} . 

-p 

Clearly, from Proposition 13. 1 1| 3|) (symmetry) it follows that = ^5^/4 • Together with (|5.4p . this 

yields that 

(5n) |^ W4 (t,t)-L((0,0),(t,t))| P ^ 

Here we use the fact that for any 77 > 0, with probability one, max |z 57I y 4 (t), Z*__^ 4 (t)| < r/t for t 

large enough. Then we can use stationarity and the fact that Ezv a (l) < +00 to conclude that for 
any e > 0, there exists 77 > 0, such that P(f 57r / 4 (77t) > ei) < e. 
By Corollary 14.21 we get 

EL, W4 (M) = E(L, W4 (t,t)-L, W4 (0,t))+EL ! , W4 (0,t) 
= 2tE(i^ r/4 (l)). 

The sub-additive ergodic theorem applied to L 57r / 4 (t,i) implies that, with probability one, 

hm = 2Ei/ 5jr/4 (l) . 

Combining this with Theorem 12.11 and (]5.1ip . one gets (|5.9p . 

Now we need to address the uniqueness of v a . Suppose v G is ergodic and time invariant. 
Define = Z_,(i, i tan(a)) and Zh(t) = Z v (t + h, ttan(a)). Now define 

Z(t) = argmax + L((—t + 2, — t tan(a)), 0)) 

z<t 
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and 

Zfi(t) = argmax (u(z) + L({—t + z, — t tan(a)), (h, 0))) . 

z<t+h 

Here, we take the right-most location of the maximum. The intuition for Z{t) and Z^if) is 
that we place the origin at (— t, — t tan(a)), and look at the exit-point for the path that starts 
at (— t, — itan(a)), picks up mass from v and then goes to 0, resp. (h, 0). Clearly, we have 

(Z(t),Z h (t))^(Z(t),Z h (t)). 

Since v is ergodic, and by our choice of a, v satisfies (|5.5p . Since the proof of Lemma I5T21 uses a 
Borel-Cantelli type argument, it is not hard to see that we can use the same ideas to prove 

(Z(t),Z h (t))/t^ (0,0). 

This means that the two paths w((—t + Z(t), — t tan(a)), 0) and w{{—t + Z^it), — t tan(a)), (h, 0)) 
will converge in any bounded box to the a-rays H7 a (0) and va a ((h, 0)) respectively (this follows 
from Theorem 12. 2( 2)). However, these two a-rays will coalesce, which means that with probability 
1, there exists to > such that for all t > to, the two converging paths coalesce, which in turn 
implies that Z{t) = Z^it) (because they are both the right-most point where the maximum takes 
place and, as soon as they coalesce, they get the same exit point). Now define 

L(t) = sup (u(z) +L((-t + z, -itan(a)),0)) 

z<t 

and 

L h (t) = sup (u(z) + L((-t + z, -ttan(a)), (/i, 0))) . 

z<t+h 

We also have that 

(L(t),L h (t)) = (L u (t,ttan(a)),L u (t + /i, ttan(a))). 
Furthermore, if t > to, then 

L h (t) - L(t) = B a (0, (0, h)) = v a ((Q, h]). 

This proves that 

M r t an(a)((*^ + /l ]) = + K 1 tan(a)) - L u (t, t tan(a)) 

A u a ((0,h]). 
Since v is time invariant and ergodic, we see that 

u((0,h]) 2 M t \ an(a) ((t,t + h}) ^ u a ((0,h]). 

In principle, we need to show convergence for a finite number of /i's simultaneously, but it is not 
hard to see that the ideas we used can be extended to that case, at the cost of some notational 
burden. Note that we have proved that for any deterministic v satisfying (|5.5p . ^(t anQ ,v([^ t + h]) 
converges in distribution to v a {h), as a process in h. This shows that in a rarefaction fan, the fluid 
process converges locally to the correct equilibrium process (local equilibrium). □ 

Corollary 5.4. 

E(*£(l)) = — jffl=. 

2\J tan(a) 
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In particular, for all a G (ir,3ir/2), we have 

%(1)).EK(1)) = ^. 

Proof: Remember that v^{ x ) = B a (Q, (0, x)). For a = 57r/4, the result follows Proposition I3.1t [3j) 
(symmetry) and Theorem 15.31 Now use the map (x,t) \— )■ (px,t/p) to see that 




□ 



For the classical Hammersley model, we know that if v\ is a Poisson counting process of intensity 
A, then v\ is time invariant and ergodic. Therefore it must be equal in distribution to v a , for some 
a G (vr, 3tt/2). We also know that in the classical Hammersley process, 

L Px (0,t) = i? 1/x ((0,t}). 



So for A = 1, we get 
Since we know that for any a G (n, 2>n/2) 



L Pl (0,t) = vi((0,t]). 



L Va (0, t) = B a (0, (0, t)) = B 5w/2 _ a (0, (t, 0)), 
using I3.1t f3j) , we can now use Corollary 15.41 to conclude that v^u = v\ ■ Consequently, 

1 = EP!(1) = Ei^ /4 (1) = ^Vtan(5vr/4) = ^ , 

which proves that 7(^1) = 2. We remark that the proof that the Poisson process is time invariant 
does not depend on the value of 7(^1). It only relies on an explicit calculation of the generator 
associated to Mj. (See Lemma 8 of Aldous & Diaconis [1], or Theorem 3.1 of Cator & Groeneboom 
0) 

Corollary 5.5. In the classical Hammersley model, we have that 7(1) = 2 and that v a = ^a(q) 
where A (a) = Vtan a. 

For general weight distributions F, we were not able to get more information on v a (not even 
a guess for a good candidate). In particular, we do not know how to calculate j(F). This does 
seem to be the most important contribution of the interacting fluid representation: once we have 
a good candidate for v a , we can check it by showing that it is invariant under the evolution of the 
interacting fluid. In fact, even in the results for the classical Hammersley case found in Aldous & 
Diaconis [Ij and Cator & Groeneboom [TJ [8] , this is where the interacting particle process proves 
its worth. 

5.1. Mixing property of v a . We will show that the measure v a has the following mixing property, 
usually called strong mixing in dynamical systems. We consider the a-algebra T = a{u a ((a, b]) : 
a < b G M} on the sample space Q, defined by the compound Poisson process. We can define the 
translation as an ^-measurable map from to 0, simply by translating all Poisson points by 
the vector (t,0). 

Proposition 5.6. For each a G (ir,3ir/2), v a satisfies 

(5.12) V4,BeJ: lim F(Ar\T7 1 (B))=F(A)F(B). 

t—too 

In particular, this implies that v a is ergodic. 
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Proof: From translation invariance and a standard approximation of sets in J 7 , it is enough to prove 
(|5.12p for all A, B € Fu ■= a{v a ((a,b}) : a < b G [0,/i]}. Consider the paths zu((-t, -t tan(a)), 0) 
and t, — t tan(a)), (h, 0)). Almost surely, these paths will converge to w a (0) and w a ((h, 0)), 
respectively, on any finite box. This means, that if we define for a,b G [0, h] 

i/W((a, b\) = L((b,0), (-t, -t tan(a))) - L((a, 0), (-t, -ttan(a))), 

then for t big enough, we have &4 = fa\[o,hy Clearly, r~~} h (B) is independent of Va\ since they 
depend on the Poisson process to the left respectively to the right of the line {— t} x R. Define the 
event 

C t = {Vs>t: v& = v a \ m } 

and denote A® the counterpart of the event A in J 7 ^ := a{i>a\(a, b]) : a < b G [0, h]}; although 
it is intuitively clear what is meant, we will make this more precise at the end of the proof. Then 

|P(Anr^(B))-P(A)P(iJ)| < \¥(Anr- + \(B)nC t )-¥(A)¥(B)\+F(Ct) 

= \F(A® n r-\(B) n C t ) - F(A)F(B)\ + P(Cf) 

< |P(AW n r f -\(B)) - P(A)P(B)| + 2P(C t c ) 
= P( J B)|P(A ( * ) ) - P(A)| + 2P(C t c ) 

< 4P(Cf). 

The proposition now follows from the fact that P(C t c ) C0. 

To see what is meant by AW, we define the index-set / = {u a ((a,b]) : a < b G [0, fo]} and 
/(*) = {v®((a,b]) : a < 6 G [0,/i]}. There is a canonical bijection «:/—)■ jw. Define ,6 as the 
product d-algebra on R J , and likewise . Extend the canonical map i such that i : R J — > R^ (t) . 
Define the map 

: O -> R / : w ^ {i/ a ((a,6])(w) : a<be[0,h]}, 

and likewise ^ : Q ->• M /(t) . We know that J), = <£ _1 (,B) and F$ = ^(B®). This means that 
there exists U G B, such that A = ^(U). We define A® = ^ (i(U)) . □ 



6. Central limit theorems for the Busemann function in the classical model 

Our geometrical approach yields a very explicit description of the fluctuations of the Busemann 
function in the classical Hammersley model. We first notice the following relations that will be 
derived from Proposition 13. II and Proposition 13.21 

Proposition 6.1. Consider the Hammersley last-passage model with random weights and recall 
(3 := (cos {5, sin a). 

• If (3 G [0,7r/2] then, as processes, 

B a (0,-P) = L„ a (-P); 

• If 13 £ [tt/2, 7t] then for each t>0 

B a (0,tf3) = v* (isin/3) - v a (-tcos/3) . 



We note that, in the classical model, we have independent increments even if the probability of the event Ct does 
not decay to very fast. This indicates that, to show mixing by using these events may not be the best strategy. 
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Proof: The first statement follows directly from Proposition 13.21 To obtain the second relation 
note that, by additivity and anti-symmetry (Proposition 13, ip . 

B a (0,t/3) = B a (0, (t cos 0,0)) +B a ((t cob p,0),tj$) 
= B a {(t cos 0,O),tff) - B«((tcoe/9,0),0) . 
Now use translation invariance, translating all four points by — (icos/3,0). □ 

In the classical set up, Corollary 15.51 together with Proposition 16.11 implies that: 

(1) if £ = then 

B a (0,tj3) = X a (tcos(3) , 

where X a is a Poisson process of intensity A (a). 

(2) if e (0,vr/2) then 

B a (0,-P) =U M , (•/?), 



where A = A(a) = \/ta,n a; 

(3) if € [it/2, it] then for each t > 

B a (0,tj3) = Y a (t sin f3) - X a (-tcosf3) , 

where Y a and X a are two independent one dimensional Poisson processes variables of in- 
tensity 1/A(a) and A(a), respectively. 

Baik and Rains (2001) proved the following central limit theorem for the Hammersley classical 
model with external "sources". Let $(x) be the standard normal distribution function, and let 
Fq{x) be the zero mean Tracy- Widon type distribution function introduced in Definition 2 of [3]. 

(4) If A G (0, 1) then 



. LpJs.s) - (1/A + X)s \ ^. . 
lim P xy \ ' K l —i- < x = $(s) : 

- 1 (7T7a^a) s v2 - ^ h 



(5) If A = 1 then 

(6) If A > 1 then 



{'■"^l 2s < x ) - /■;,(,,-,: 



s->oo 



- 1 (7a^T7a),V2 - I ^ 



S— i-OO 



These results naturally lead us to a central limit theorem for the Busemann function 

B((3,-) := ^4 (0, (• cos /3, • sin /?)). 

By (|2.4p . w.l.o.g. we can restrict ourselves to a = 57r/4. In this case we have 

RB(p,t) = (cos/3 + sin/3)t. 

Corollary 6.2. Consider the Hammersley classical last-passage model. Then 

( B[tt/A,£) -y/2t \ 
fe P 2-1/4/3 <* =*b(*), 
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while for j3 G [0,7r/4) 



and for /3 G (vr/4, ir] 



, B(B,t) - (cos/3 + sin/3)i \ 
lim P y ,\ > V — , 1/0 < x = $ x , 



t— >oo 



(cos /3 + sin #)t \ 
1 (v , sin ( 0-cos/^)t 1 /2 - i v ' 



t— >oo 



Proof: The statement for /3 = and /3 G (7r/2, 7r] follows (1) and (3), together with the central 
limit theorem. Now, let us take (3 G (0, n/2). Consider the map (|2.4p with p{j3) = -^/tan /3. Then 

B((3,t) = Lp p (s,s), where s = (\/s'm /3 cos and for /3 G (0,7r/4] we have p G (0,1] while for 
(3 G (vr/4, 7r] we have p G (l,oo). □ 



6.1. The crossing formula. In the classical model, the exit point formula for the equilibrium 
regime, proved by Cator and Grooeneboom [8], is 

(6.1) VarLp A (x,t) = -Ax + ^ + 2XEZ Px (x, t)+ , 

A 

where VarX is the variance of X, and X + := m&x{X, 0}. Together with (2) and Proposition 15.11 
this relates the variance of the Busemann function at (x, t) with the position of the crossing point 
of the a-ray starting at (x,t): 

(6.2) Var5 Q (0, (x, t)) = -(\Ztana)s + . = + 2(Vtan a)EZ a (x, t)+ . 

Vtana 

In particular, 

Var5 57r/4 (0, (t,t)) = 2EZ &r/4 (i,i)+ . 

We note that this crossing point formula can be seen as a version of the scaling identity £ = 2%, 
where x and £ are the critical exponents that measure the order of magnitude of the fluctuations 
of Busemann functions and crossing points, respectively. 

7. The multi-class process and second class particles 

For two positive measures v and v on M, we say that v dominates u, notation u > u, whenever 
> V {I) f° r a ll measurable / CI, 

Proposition 7.1. Suppose we have two measures v, v G M such that v > v. Define the corre- 
sponding interacting fluid system as and M" , using the same weighted Poisson process (basic 
coupling). Then Mf > M\ (as measures). If 'P(z) = u(z) for all z < 0, then Mf ([0, x]) - M?([0, x]) 
is non-increasing in t for all x > 0. 

Proof: Fix an interval [-K, K] and a time t. There exists (a random) M > such that 
and Mj* restricted to [-K, K] only depend on Poisson points in [— M, K] x [0, t] and on v and v 
restricted to [— M, K] (it is not hard to see that we can take M = Z u (—K,t)). This means that 
we are only dealing with a finite number of Poisson points, so if we can prove that the premise 
li Mg > Mg for all s < t" implies that > M", we will have proved the first statement, since it 
is obviously true for t = 0. Suppose there exists a Poisson point at (xo,t) with weight oj for some 
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xq G [—M,K], since otherwise the implication is immediate. We then know, using Proposition 14.21 
and (|4.ip . that if xo < x < y, 

M?((x,y]) = (M^_((x ,y])-u) + -(M^((x ,x])-u) + 
< (M?_((x ,y])-u>)+ - (M£((x ,x]) 
= M?((x,y]). 

The inequality follows from the fact that if A > B and A > B > 0, then (A — u>)+ — (B — w) + < 
(A + A — w) + — {B + B — If x < xo < y or x < y < xq, the implication is straightforward, 
following a similar split up. 

The second statement follows from a similar reasoning: suppose there is a Poisson point at (xo, t) 
with weight to. If x > xo > 0, 

Mf((x ,x])-Mn(x ,x]) = (Mf_((x ,x])-o;) + -(Mr_((xo,x])-a;) + 

< Mf„((xo,x])-M t ^((x ,x]). 

The inequality follows from the fact that [A — c)+ — [B — c)+ < ^4+ — B + whenever c > and 
^ > B. Since Mf ([0,x ]) - M^([0,x ]) = Mf_([0,x ]) - M^_([0,x ]), this shows that 

Mf([0,x]) -M t "([0,x]) < Mt([0,x])-M t 1([0,x]). 

Now suppose xo < and x > 0. Note that under the condition on u, we have that for all s > 
and all e > 0, Lp(—e, s) = L u (—e, s), so 

Mf ([0, x]) - M;([0, x]) = L p (x, s) - L„(x, s). 

When Lp(x,t) does not use the weight at (xo,i), we know that Lp(x,t) = Lp(x,t—) and that 
L u (x,t) > L u (x,t—), which implies the desired result. If Lp(x,t) does use the weight at (xo,t), 
then it is not hard to see that L u (x,t) will also use the weight at (xo,t) (the longest path corre- 
sponding to v is always to the right of the path corresponding to v), which means that only the 
mass on the x-axis strictly to the left of is used, and therefore M t "([0, x]) = M^([0,x]). Finally, 
when x = x or x > x, we get that Mf_([0, x}) = Mf ([0, x]) and Mf_([0, x}) = Mf([0, xj). □ 

In other words, Proposition 17.11 tells us that the interacting fluid system is monotone: if one 
starts the fluid process with the same Poisson weights (basic coupling) and with ordered initial 
configurations, then the order is preserved for all t > 0. This coupled process is called the multi- 
class fluid system. The multi-class system is just a convention to describe a coupled process with 
ordered initial configurations (Ferrari & Martin [10j). 

7.1. The multi-class invariant process. With Theorem 12.21 in hands, for any countable D C 
(-7T, 3-7r/2), one can construct simultaneously a collection of equilibrium processes \v a : a £ D} by 
using the same Poisson weights on R x R_ and the Busemann functions B a . It turns out that this 
collection respects the order induced by the angles a £ D. More precisely: 

Theorem 7.2. If a > a then v & > v a . In particular, for any countable subset {oti : i G Z} C D, 
if one runs simultaneously (basic coupling) the interacting fluid processes on R x R + with initial 
measures (y ai : i G Z) then, whenever cti > ctj, M^ ai > M t " J for all t>0, and 

{M" ai : % G Z) = (u ai : i G Z) . 
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Proof: Let z' > z £ R. Let m be the crossing point between tu s ((z,0)) and ti7 a ((V, 0)). Further- 
more, denote c as the coalescence point of the two a-rays w a ((z,0)) and zu a ((z', 0)), and denote c 
as the coalescence point of the two a-rays cc7q((z,0)) and ^((z', 0)). Then 



'a(M)-fa(M]) = {L(c,(z',0))-L(c,(z,0))} 

L(c,(z',0))-L(c, 0,0))} 



= L(c, 0)) - {L(c, m) + L(m, 0))} 

+ L(c,(z,0)) - {L(c,m) +L(m, (z,0))} > 0. 

Notice that the Busemann functions B a are a function of the compound Poisson process (P,w): 
-Ba(') •) = B a (-, -)(P,w) (here w denotes the weights). Since p + P, the translated version of P, 
has the same distribution as P, we get that 

{B ai (-,-)(P,w) : i€Z}i{S a .(-,.)(p + P,w) : ieZ}. 
This shows time invariance for the Busemann multi-class process. 

□ 

This result is also new in the classical Hammersley interacting system, where a different and 
explicit description of the invariant process with a finite number of classes is given in Ferrari & 
Martin [TO] . 

7.2. Law of large numbers for second-class particles. Proposition 17.11 can be used to define 
the notion of second-class particles. In the interacting fluid system we can define it analogously to 
the interacting particle case, with a slight adaptation due to the continuous weights. We start by 
changing v into u, by putting an extra weight e > in 0, so 

v([0,x]) = v([0,x]) + e for x>0. 

With this new process, and using the same Poisson weights, we define M^. Clearly, 

M*([0,x])<M?([0,x])+e. 

Now define the location of the second class particle X v (t) as 

X u {t) = inf{x > : Mf(x) = M v t (x) + s} . 

By Proposition 17. 1\ X u (t) is a non-decreasing function of t, meaning that the second class particle 
moves to the right. In fact, the extra mass e will spread out, and our definition coincides with the 
rightmost point of this spread-out mass. This is a natural choice, since we will show that it does 
not depend on the total mass e, while for example the leftmost point does depend on e. 

There is the following important connection between the longest path description and the second 
class particle. Let v + be the process defined by v + (x) = v{x) for x > 0, and by v + {x) = — oo for 
x < 0. We also define the process v~ by setting v~(x) = for x > 0, and v~(x) = v{x) for x < 0. 
Then 

r r a.\ f sup{L((z,0),(x,t)) +v(z) : 0<z<x} ifx>0 
[ — oo if x < 0, 

and 

L u - {x, t) = sup{L((z, 0), (x, t)) + v(z) : z < and z < x}. 

Clearly, 

L„(x,t) = max{L u +(x,t),L u ~(x,t)} . 
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Now suppose x > 0. If L v +{x,t) > L v -(x,t), there exists a longest path that does not use any 
weight of v on (— oo,0). This means that if we add a weight e > in the origin, Lp(x,t) = 
L p +(x,t) = L v (x,t) + e. Using Proposition 14.11 we see that this means that M^{x) = M"(x) + e, 
so X u (t) < x. If on the other hand we start with X u (i) < x, we conclude that M"(x) = Mf(x) + e, 
using Proposition 17.11 and the fact that M" and M" are right-continuous. This in turn means that 
Lp(x,t) = L u (x,t) + e, which is only possible if L u +(x,t) > L v -{x,t). We have shown that 

(7.1) {X u (t) <x} = {L u + (x, t) > L„- {x, t)} . 
Note that this can be rewritten as 

(7.2) {X v (t)<x} = {Z v (x,t)>0}. 

This means that the path of the second class particle corresponds to a competition interface, a fact 
well known for the totally asymmetric exclusion process (Ferrari & Pimentel 2005). This allows us 
to show that the second class particle satisfies a strong law whenever v + and v~ have asymptotic 
intensities. The proof of this does not use a coupling of two invariant versions of the fluid process, 
as is usual in the interacting particle case, but it uses the longest path description in a direct way. 
We would like to point out that in our general set-up, with random weights on the Poisson points, 
we do not have an equivalent of Burke's Theorem. This means that the time-reversed process is 
not a Hammersley interacting fluid system. Therefore, the path of a second class particle in general 
does not coincide in law with a longest path in the interacting fluid system, in contrast to the 
classical case, where the statement is true. However, we do have the following connection. 

Proposition 7.3. Assume that the distribution of v is translation invariant. Then, for any t > 0, 
we have that 

X u (t) - x = -Z l/ (x,t). 

Proof: This follows almost immediately from (j7.2|) . since that equality can be rewritten as 

{X u (t) - x<h} = {Z u (x + h,t) > 0}. 
Now use translation invariance to see that 

Z v (x + h,t) = Z u (x, t) + h. 
Combining these two equations proves the proposition. 

□ 



When we consider all a-rays starting at the line R x {t} and we move from left to right, (X Ua (t),t) 
is the first point where the a-ray passes the origin. It is tempting to think that the a-ray starting 
at (X Va (t),t) actually passes through the origin, but this is false in general. In fact, after time t, 
most a-rays will have coalesced with other rays, and the crossings with the x-axis will be quite far 
apart; we would conjecture they are order i 2 / 3 apart. 

Proposition 17. 31 allows us to use Theorem 12.11 and Lemma [5 .21 to prove a strong law for the second 
class particle in the case of v a . However, we are able to prove a strong law even for deterministic 
initial conditions that satisfy a density property: 

Theorem 7.4. Assume that 

(7.3) lim - - = lim - - = — Vtana . 

i->oo x y— oo X 2 

Then, with probability one, 

X u (t) 1 



lim 

t^oo t tan a 
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Proof: It is enough to prove the statement for a = 5ir/A. Suppose e > and X u (t) < t — let. 
Define n = \t\. Then for t large enough, we have 

X u {n) < (1 - 2e)n + 1 + 2e < (1 - e)n . 

By (j7.2[) . this implies that Z u ((l — e)n, n) > 0. In the same way as in Lemma 15.51 this can happen 
only for finitely many n > 1, which gives 

liminf^l > i. 

t— >oo t 

Bounding the limit from above can be done using the analogous argument. 

□ 



8. Lattice last-passage percolation 

In the lattice last-passage percolation model we have a collection {w z : z G Z 2 } of i.i.d. non 
negative random variables indexed by lattice points z G Z 2 . In this set up, one can define last- 
passage times for x < y (we put an i on the top to indicate that it refers to the lattice model) by 
maximizing over up-right paths connecting x to y: 

L e (-x-,y) ■= max {V]w x /}. 

roen x,y) 

Similary to the preivious case, one can start with a non-decreasing process z/ = (u (z) : z G Z) 
and define the evolution by 

M?{{x,y\) := L u e(y,t) - L u ,(x,t) , (x,yGZ, t > 0) 

where 

Zv(x, t) := sup \v e (z) + L e ((z, 0), (x, t))\ 

(i/ must have a left density analogous to (|4.2p ). For a survey in lattice last-passage percolation 
(and its connection with the totally asymmetric exclusion process) we address to |18j . 

If oj z has an exponential distribution of parameter one, then the limit shape is given by (y/x+s/i) 2 
and a result similar to Theorem 12.21 holds (Pimentel & Ferrari [12]). This allows us to construct 
Busemann functions B l for a G (tt,3tt/2). The same method developed to prove Theorem 15.31 can 
be applied to this case. Since for exponential weights we also know the invariant measure, we have 
that: 

Theorem 8.1. For a G (7r,3vr/2) let 

u e a (x) = B e a ((0,0), (x,0)) for xGZ. 
Then i/~ is the unique ergodic process on Z that satisfies 

Mf a = i/ a for all t > . 

In particular, for any x,y G Z, M t a (x,y) is distributed like a sum of i.i.d. exponential random 
variables of intensity 

\/cos a 



U) p(a) 



^cos a + V sin a 
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We remark that p(a) is obtained by solving 

/H / ,o 1 tana 

(1 + Vtan a) = — H . 

P 1 -P 

The left-hand side corresponds to the limit shape of the lattice model, in direction (1, tan a), while 
the right-hand side corresponds to summing up, in the same direction, the expectation of the 
equilibrium process in the horizontal and vertical axis. (Compare this with (|5.9p and its proof.) 
The dual process on the vertical axis, denoted by z/^*, is distributed like a sum of i.i.d. exponential 
weights of intensity 1 — p(a); see also Balazs, Cator and Seppalainen [4]. 

To develop our method in the model with general i.i.d. weights, one would need an inequality 
similar to (|5.6p for the respective limit shape. This is, however, one of the most challenging problems 
in lattice last (and first) passage percolation models. 

9. Appendix 

9.1. Proof of Proposition 13.11 The translation invariance of the underlying compound two- 
dimensional Poisson process, 

K:qe(P + z)} = {w p :peP}, 

implies (P). 

Anti-symmetry follows directly from the definition of the Busemann function. Now, by taking a 
coalescence point c between w a (x), w a {y) and w a (z) we have that 

B a (x, z) = L(c, z) — L(c, x) 

= L(c,z) - L(c,y) + L(c,y) - L(c,x) 
= B a (y,z) +£ a (x,y), 

which clearly shows additivity, and finishes the proof of ([2]) . 

Define S as the reflection in the diagonal x = t. Then S(P) has the same distribution as P, and 
x < y 44> S(x) < S(y). This shows that for all x,y£ R 2 , 

L(xj)^(5(x),S(y)). 

This clearly proves ([3]). 

Since L(c,y) > L(c,x) whenever c < x < y we have that B a (y,x) > whenever x < y. Note 
that 

B a (x,x + (1,1)) > L(x,x + (1,1)). 
Clearly, EL(x, x + (1, 1)) > 0. Also, 

E J B Q (x, x + (1, 1)) = E J B a (x, x + (1,0)) + E J B a (x, x + (0, 1)). 

Now the first three properties prove that MB a (x, y) > 0. 

To prove that its expected value is finite, without loss of generality, assume that y = 0, that 
a = 57r/4, and that x = (—1, —1)- Let z be the intersection point between ro 57r / 4 (0) and the one 
dimensional boundary of { p £ I 2 : p<x}. By taking a coalescence point c for and x such 
that c < z, we have that 

£ 57r/4 (x,0) = L(c,0) -L(c,x) 

= L(c,z) +L(z,0) - L(c,x) 
= L(z,0)-{L(c,x)-L(c,z)} 
< L(z,0). 
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Now assume that r\ > is a small constant and that |z| < r]r. Then 

L(z, 0) < max{L ((-1, -rjr), 0) , L ((-ryr, -1), 0)} < P(rjr) 

where P (r) denotes the number of Poisson points in the strip composed of points (x, i) such that 
x G [-1,0] and t G [-r,0], or i G [-1,0] and x G [-r, 0]. Hence 

P (L(z, 0) > r) < P (|z| > r?r) + P (P(r/r) > r) . 

Denote by Co((— 1, — 1), 9) the cone with axis the half-line starting at through (—1,-1), and 
angle 29. Lemma 2.3 of [9] (5-straightness of geodesies), states that for S G (0, 1/4), 9 = tt/5 and 
M > big enough, there exist Co, c\ > such that for all p G Co((— 1, —1), 9) and |p| > M, 

P({q£i 2 : p G w(q,0)} G: ^(p^lpl- 5 )) < c e- Cl l p l K . 

If \z\ > r]r, then there must be p G 3Co((— 1, — 1), 9) with |p| > rjr, such that p G vj^/^O). 
However, for r big enough, w^/^O) will leave the cone Co(p, c|p>| — 5 ) , since it must follow the 
(negative) diagonal, implying 

{qel 2 : p G ro(q,0)} £ Co(p,c|p|^). 

Therefore, P(|z| > r/r) is integrable in r (for all rf). On the other hand, by choosing rj small enough, 
one can make P (P(r/r) > r) integrable over r > as well. This finally implies that EL(z, 0) < oo, 
and consequently, E£> 57r/ / 4 (x, 0) < oo. 

By additivity and anti-symmetry, to prove cadlag, we can restrict our attention to x, y, p = 0, 
A ^ and q = (1,0) so that we are varying in the horizontal direction close to the origin. For 
q = (0,1) the argument is similar. Fix eo > and assume that B a (0, (eo, 0)) > 5. Then one 
can find Co such that the path w(cq, (eo,0)) is disjoint from the path w(cq, (0,0)) and contains at 
least one Poisson point p° = (p?,^) with p\ G (0,eo). Define the set A as the region enclosed 
by tu(co,0), the x-axis and the y-axis. Now take e\ = Pi/2 and assume that B a (0, (ei,0)) > 5. 
Then, similarly, one can find ci > Co such that the path ro(ci, (ei, 0)) is disjoint from the path 
vj(ci, (0,0)) and contains at least one Poisson point p 1 = ip\,p\) with p\ G (0,ei). But, now, we 
must have that p 1 G A. In this way, if liminf e _ i .o B a (0, (e, 0)) > 5, then one could find infinitely 
many Poisson points within the bounded set A, which contradicts the fact that, a.s., the Poisson 
process is locally finite. □ 
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